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Executive Summary 


This report summarizes the efforts performed under AFOSR contracts FA9550-09-C-02I3 and 
FA9550-O7-C-O039, which constitute one project. The project includes the two major 
components: modeling of the surface dielectric barrier discharge (SDBD) and CFD studies 
SDBD flow control capabilities. The first component has been focused on developments of self- 
eonsistent theoretical models and robust computational tools providing adequate from-first- 
prineiples simulations of heat and momentum sources induced by surface discharges. The 
physical model and numerical code, which were developed to predict SDBD process in 
atmospheric air, allow us to simulate the total SDBD cycle including the discharge formation 
phase and the discharge relaxation phase. The predieted values of the discharge length and the 
surface charge density agree well with the available experimental data for positive and negative 
eleetrode polarity. This validates the developed physical model and computational code. 
Numerical solutions clarified basie features of the discharge evolution in space and time. 

Using the numerical solutions and available experimental data we have developed a physics- 
based phenomenological model, which sheds light on physical processes governing the ignition 
and eeasing of mierodiseharges for the negative going half-cycle of applied alternating voltage. 
The obtained analytical formulae capture dependencies of the body foree and discharge length on 
the SDBD parameters. Using the body-force values predicted by the phenomenological model 
and numerical results for the spatial distributions of momentum and heat sources, we have 
obtained analytical approximations of the SDBD-indueed body force and heat source. The source 
terms were incorporated into Navier-Stokes solver. This computational tool has been used in the 
seeond component of the project dealing with SDBD applications to flow control. 

It was shown that SDBD actuator can strongly affect the laminar boundary layer on a flat plate. 
Depending on distributions of applied voltage it is feasible to generate tangential jets or 
concentrated vortices in the near-wall flow. This flexibility and very short characteristic 
timescales make SDBD actuators attractive for control of viscous flows associated with the 
boundary-layer separation, laminar-turbulent transition or a combination of both. 

CFD studies of the vortex flow past a delta wing at high angles of attack have been carried out to 
estimate feasibility of the vortex flow control using SDBD actuators. Computations were 
performed for a delta wing with sharp leading edges of 60-degree sweep angle at the free-stream 
Maeh number 1.5. It was shown that, for this configuration, the flow control strategy should be 
foeused on the vortex breakdown. Computations have been performed for the wing-apex SDBD, 
leading-edge SDBD and multi-element SDBD actuators. The heat and momentum sources 
produced by these actuators were symmetric with respect to the wing centerline. It was found 
that the vortex-burst loeus can be controlled by the aforementioned actuators. Namely, the 
breakdown point moves downstream with the SDBD acting in the downstream direction while it 
moves upstream with the SDBD acting upstream. However, the integral aerodynamic forces are 
weakly affected. The actuator causes about 3% variations of the lift and drag coefficients. It was 
also shown that an asymmetric forcing, which can be produced by SDBD actuator on one side of 
the wing surface, triggers strong oscillations of the vortex burst loci. I his. in turn, causes 
pulsations of the rolling moment. Because of this detrimental effeet it is not clear if SDBD 
actuators are robust for the vortex flow control on supersonic delta wings at high angles of 
attack. Further numerical and experimental studies are needed to clarify this issue 
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1. Background 

Our original effort (2001-2003) was devoted to theoretical studies of surface discharge plasmas 
and other plasma devices to quench stall-slip departure due to asymmetric vortex shedding on 
aircraft nose tips. The application is to augment directional control/authority, agility and provide 
tailless capability of high-speed flight vehicles [I]. We started with theoretical modeling of 
vortex symmetry breaking on a sharp circular cone at high angle of attack From slender body 
theory, we conceptualized a theoretical model of the vortex dynamics and estimated feasibility of 
vortex control by surface discharges [2]. Within the framework of an idealization of the large- 
scale vortex dynamics as an inviscid flow consisting of a pair of straight feeding sheets 
culminating in point vortices that idealize the feeding sheet rollup, we found that the symmetry 
breaking phenomenon is associated with nonlinear instability of a vortical singularity saddle 
point occurring as it rises above the body surface. Results of this effort are summarized in [2,3]. 

These theoretical studies were complemented by experimental effort [4,5] of the Institute of 
Theoretical and Applied mechanics (1TAM, Novosibirsk, Russia) that was performed under the 
ISTC project 2235. The experiments were conducted in ITAIVTs T-324 subsonic w ind tunnel on 
a slender cone equipped with spark discharge devices. These devices created localized plasma 
heating that strongly displaced the cone boundary layer, triggering early separation. In accord 
with the theoretical predictions of [2,3], these effects stabilize the symmetric flow pattern and 
avoid symmetry breaking. The experiments [4,5] confirmed our concept of plasma stall-spin 
departure control in which primarily inviscid large-scale vortex dynamics can be controlled by 
spark discharge relocation of the boundary layer separation loci that are the origins of the feeding 
sheets. Our electromechanical system that implements this principle of viseous-inviscid 
boundary layer interaction control has been patented by Rockwell Scientific (RSC)-Russian 
team 3 Because of the importance of the application of this technology and the need for further 
physics-based modeling, we have continued the previous effort with emphasis on plasma How 
control actuators such as spark discharges (SO) as well as dielectric barrier discharges (DBD) 
that fall into a wider class of discharges known as surface barrier discharges (SBD). 

In 2004-2006 wc have improved the fidelity of our simulation of the viscous-inviscid interaction 
between the large-scale vortex dynamics and boundary layer that leads to symmetry breaking. In 
connection with this, a conical Navier-Stokes (CNS) approximation was applied to model vortex 
dynamics over a sharp cone at high angles of attack [6]. This simulation captured the symmetry 
breaking phenomenon. Special attention was paid to modeling of the role of turbulence in the 
symmetry breaking. It was shown that our CNS solver reproduces symmetric and asymmetric 
vortex fields and can be used for predictions of the critical angle of attack at which the vortex 
symmetry destabilizes. Using this tool we carried out parametric studies of surface plasma 
discharge effects on the vortex structure [7] It was shown that CNS predictions are qualitatively 
similar to the experimental observations [4]. 

In connection with rapidly growing technologies relevant to flow control of high altitude flight 
vehicles. Dr. Norman Malmuth suggested to focuse our 2007-2010 effort on plasma How control 
for supersonic delta wings at high angles of attack. Similar to conical forebodies, the boundary 
layer globally separates from the leeward surface of delta wing at high angles of attack (a) [8]. 


'Malmuth, N., Fedorov, A., Soloviev, V., Maslov, A., Zharov, V., Shalaev, I., “ Surface Plasma 
Discharge For Controlling Forcbody Vortex Asymmetry,” US Pat. No. 6,796,532 (09/28/2004) 
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The flow field contains two strong vortices generated by the rollup of the shear layer emanating 
from the separation lines near the wing leading edges. Large radial and axial velocities of these 
vortices reduce the pressure on the wing surface under the vortex axis. This leads to increasing 
lift. With growing a, the lift increases nonlinearly. For sub-eritieal angles of attack. a<a cr , the 
vortices are globally stable. This makes the flow field symmetrical and only weakly sensitive to 
flow perturbations. For supercritical angles, a>a cr , the vortiees exhibit a large-seale instability 
that is characterized by rapid deceleration and eventual stagnation of the vortex core flow as well 
as growth of the vortex eore eross-seetion These events are accompanied by strong oscillations 
and eventual breakdown of the vortieal structure [8-10]. This phenomenon, ealled vortex 
breakdown or vortex burst, leads to a sudden decrease of the lift-to-drag ratio (L/D). If the vortex 
breakdown on one wing side oeeurs earlier than on the other, then yawing and rolling moments 
are induced due to flow asymmetry. 

It is well known that leading-edge vortex flow is extremely sensitive to variations of surface 
shape and/or external forcing near the separation line. For wings with sharp leading edges, the 
separation line practically coincides with the leading edge. For rounded leading edges, the 
separation line has a complex curved shape and runs along the leading edge from the apex to the 
wing tip. Changing the surface boundary condition in quasi-steady state (for example, expanding 
surfaee curvature) near the separation line ean eause a global ehange of the developing vortex 
[9]. Loeal instabilities can be also exploited for manipulations of separated vortices. As shown in 
[11] very small foreing can eontrol vortex structures in the mixing layer if the forcing frequency 
lies within the sub-harmonic frequency band. Sueh foreing ean be produced by plasma actuators 
installed near the wing leading edges (close to natural separation lines). 

Plasma discharge actuators may be advantageous compared to synthetic jets, MEMS and other 
aerodynamic schemes. This explains the rapid growth of plasma flow eontrol technology [12]. 
For this purpose, the most attractive is a surfaee dielectric barrier discharge (SI)BD) that is of 
simple construction, does not ehange its aerodynamic shape or influence the wing when it is not 
in use, allows for smooth forcing frequency and power variations, and can be used for closed- 
loop feedback control. This motivated Corke et al. [12-14] to utilize SDBD actuators for control 
of boundary-layer separation near the leading edge of a subsonic airfoil at high angles of attack. 
The experiments [15] on a subsonic rectangular wing at the chord Reynolds numbers from of 
0.35xl0 6 to ().875xl(/ demonstrated that SDBD aetuators of low power (less than 250 Watts) 
effectively suppress leading edge separation and improve wing performance at high angles of 
attaek. 

The major gap in plasma flow eontrol technology is the lack of self-eonsistent theoretical models 
and robust computational tools that ean provide adequate simulation of heat and momentum 
sources induced by SDBD aetuators. In contrast to numerous papers dealing with the volumetric 
arrangement of SDBDs, optimal surfaee configurations have been analyzed numerically only in a 
few papers such as [16,17]. Gibalov and Pietsehet [16] calculated SDBD development in air for 
constant negative polarity electrodes. In this numerical simulation, the secondary emission of 
electrons from the eathode is treated as a dominant SBD mechanism, although there are 
experimental data showing that sueh discharges ean be produced with no secondary emission 
[18]. Another shortcoming of this work (even for constant electrode polarity) is that the 
calculations were carried out for a relatively short time that did not cover the discharge decay 
phase. Accordingly, it is not feasible to compare the predicted length of the discharge region 
with available experimental data and validate numerical results. 
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The numerical simulations [17,18] addressed the more praetieal ease of a radio frequency (RF) 
surface barrier discharge. However, the modeling in [17] was performed for a diffusive discharge 
in helium at 300 Torr pressure when the discharge bums continuously and docs not have a 
streamer phase, relevant to burning in air. The work in [18] simulating SBD evolution in air with 
an alternating voltage source has the same deficiency. The authors computed all the discharge 
phases with uniform time stepping that is too large to deseribc streamer development. Without 
the streamer phase, the solely diffusive model prediets gas heating and acceleration that could be 
physieally unrealistic. Additionally, they did not account for air photoionization by UV radiation 
from the discharge region. This radiation is neeessary for streamer production. 

To bridge these gaps we have conducted combined theoretical and computation studies of SDBD 
physics. This allows us to identify key physical mechanisms associated with different phases of 
SDBD eyele and evaluate their basic parameters. Using this knowledge we have developed a 
self-eonsistent phenomenological model, whieh provides SDBD-indueed momentum and heat 
sourees for flow-control modeling. 


2 . Major accomplishments 


An overview of the projeet effort is shown in Fig. 1. The projeet includes the two major 
components: 1) plasma discharge modeling and 2) fluid dynamics modeling. 





Rtonfor hnriy Rerlangidar Vrfng Delta vuinq 

Figure 1 Project overview with major components. 
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1 he first component is focused on developments of self-consistent theoretical models and robust 
computational tools providing adequate from-flrst-prineiples simulation of heat and momentum 
sources induced by surface discharges. The second component deals with fluid dynamics 
modeling relevant to plasma flow control applications. These two components are strongly 
coupled and aimed at the development of a combined theoretical and computational toolbox for 
plasma flow control applications. 


2.1 Modeling of surface dielectric barrier discharge 


2.1.1 Accomplishments of the 2007 - 2008 effort 


Schematics of SDBD electrode layout is shown in Fig. 2. The exposed electrode is assumed to be 
loaded by a high voltage, whereas an electrode buried under the dielectric is grounded. If the 
applied voltage is greater than a threshold value, the discharge starts glowing near the exposed 
electrode edge. According to the experimental data [16], if the exposed electrode is a cathode. 


then discharge takes a diffusive form. For 
positive exposed electrode polarity, the 
discharge evolves as a streamer. The 
discharge region length for the streamer 
stage is few times greater than that for the 
diffusive stage and in atmospheric air 
approximately equals to 10 mm for the near 
threshold voltage value V % i - 5 kV, the 
dielectric thickness d zz 1 mm, and the 
relative dielectric permittivity e « 5 - 10. 



In the case of applied voltage of low 

frequency (less than 10 - 100 kHz), SDBD Figure 2 SDBD electrode layout and used coordinate 
reveals two phases of its evolution: the system, 
discharge formation phase, when plasma 

layer is created due to intense air ionization, and the relaxation phase, when this plasma decays, 
fhe discharge formation phase lasts approximately 20 50 ns. It is accompanied by a high 

amplitude electrical current pulse in the external circuit. During this phase the dielectric surface 
is seeded by charged particles (electrons or positive ions depending on applied voltage sign), 
which screen the external electric field. This stops further discharge formation and its 
propagation along the dielectric surface. The subsequent relaxation phase lasts few 
microseconds, until the electric field restores its breakdown value due to both surface charge and 
plasma decay. Accordingly, for low frequency voltage, SDBD is arranged as a set of short (about 
tens of nanoseconds) current pulses separated by relatively long (few microseconds) relaxation 
pauses. 


2. / la Statement of problem and governing equations 

Discharge simulations have been performed in 2D approximation for the electrode layout and the 
coordinate system shown in Fig. 1 Hereafter the exposed semi-infinite electrode has the height 
h e = 0.1 mm above the dielectric surface, the dielectric thickness is d 1 mm, and the relative 
dielectric permittivity is e = 8. 
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The model accounts for three types of charged particles: electrons, negative ions and positive 
ions. Their concentrations are n f , n , and n t , respectively. The analysis of kinetic processes 
[19,20] showed that negative ions are primarily represented by O ions, and positive ions by 
O^O; ions coming after a fast conversion of N 4 \N, ions in charge-exchange collisions with 
O, molecules. Positive ions 0 4 \0* arc close to equilibrium resulting only in one type of 
effective positive ion. The used simplified air plasma kinetic model was chosen to capture the 
main features of the discharge evolution, and it does not pretend to describe details of plasma 
composition. 

The transport equations for charged particles are 


dr). 7 

-h div,I 

Of 

= k t Nn ,, k r n n, k t!r n t ,n t f S fth , 

.7, = n K E. 

(2.1) 

On . 7 

-h dtvJ 

Of 

— 0.22k 1tl Nn, k df n N k t n n r 

.7 = n K E. 

(2.2) 

+divl 

dt 

= s, + S ph .S t = k\Nn, - k dr n,n, 

-0.22 k nl Nn, N . 

(2.3) 


J, = -V(D r n ( )-n t K ( E, 


where J = flux of particles, D = diffusion coefficient, K - mobility, the subscripts e.T. 

denote electrons, positive and negative ions, E = electric field, N ~ air density, 5 = electron 
souree/sink temi due to ionization, recombination and attachment processes, S - elcctron-ion 

source term due to air ionization by radiation from the discharge region The factor 0.22 denotes 
an oxygen fraction in the air, A„ kj r , Ay, k i1f and kj, are rate constants for ionization, dissociative 
electron-ion recombination, ion-ion recombination, dissociative electron attachment and electron 
detachment, respectively. 

I he transport equations arc complemented by Poisson equation for the £-field potential y 


Ay — l7re(n n n t ) for gas region ( v > 0), 
Ay 0 for dielectric layer (-d <y < 0), 

E - Vy. 

Conditions on the gas-dielectric boundary v _ 0, .v > 0 are 


0 y 

()y 

r)y 

- - 

Ox 

Ox 

n dy 

(i i 0y 


(2.4) 

(2.5) 

( 2 . 6 ) 


(2.7) 


where a = a + a t = surface charge density. Boundary conditions for the electric potential are 


y 0 at v - iL y — V for 0 < v < h C9 x < 0, 


( 2 . 8 ) 


1 1 


- arctg 

— 

2 7T 



for v —> oo, .v-> ± oo , 


(2.9) 
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where V = applied voltage (exposed electrode potential). The condition (2.9) results from an 
asymptotic solution of the equation — 0 in the semi-plane v > 0 with the boundary 
conditions ( 2 . 8 ) for d 0 . 

The rate constants are specified in Table 1, where y = redueed electric field E/N in units of 10 ’ 16 
V ein". The electron temperature T e and the ion temperature 7} are measured in Kelvin, and T, is 
calculated using the Wannier expression 

7/1 -L 

T t — T fl + \ ( 2 . 10 ) 

where T n = gas temperature, m r rn n ~ mass of ion and X molecule, k n = Boltzmann constant. 
Table 1 Ionization and sink rate constants 


Reaction 

Rate constant (cm s 

) Reference 

+ N,(OJ - 

* 2 p 4- 

N ((). ) 

A = (0.7668 + 0.018 

7)10 

r, [21], [23] 

O' 

-0 + 0 

Kx =2x10 7 

300 

[24] 




T t (K) t 


r f O; 

* 0, 

1 0, 

A;,, 1.1x10'* 

300 

124] 




J.UO 


+ 0 , - 

0 

-f 0 

- 10 . 

Kn 10 


[21] 

O + X, 

—* r 4 

X,0 

K »-2 • 

10 ‘ 

[24] 

' + (M) — A 


M (M 0 : . N : ) 

A,= 2»10' 6 (300/ T ,)' 5 

[25] 


The rate constants for the processes depending on the reduced electric field E/N together with 
A' t ,, D ,,, and T e are defined by the electron energy distribution function. They are taken from Ref. 
[22]. The analytical approximations of these rate constants are given in Ref. [21]. The ionization 
rate constants and the electron drift velocity, or the electron mobility A' t ., are further eorreeted to 
fit the data [23] (this has been discussed in Refs. [19, 20]). Finally, the following expressions for 
K e , 7V, and D t . have been used 


V s 


8GG 


7 < I 


883 (1 -I 1.26277) 1 < 7 < 20 , 


T (K) - 
D = ^ 


275-(1+ 9.G/->) 

T -+ 8045 • 7 " ,41 "'' -> < 1 

7’ -f 8645 • V’' > 1 

k„TK 


7 > 20 


( 2 . 11 ) 


( 2 . 12 ) 

(2.13) 


S 


















The mobility of (), and 0 . ions in atmospheric air is known to be 2.4 and 2.05 cm 2 /Vs, 
respectively. We imply K = 2.1 cirf/Vs for effective positive ions and K 3.2 em 2 /Vs for 
negative O ions. The electron-ion recombination and electron detachment rate constants are 
taken from Ref. [24], and the ion-ion recombination rate constant is taken from Ref. [25]. 

In the right-hand sides of Eqs. ( 2 . 1 ) and (2.3), the source term S p h is due to the photo-ionization 
of O: molecules by UV radiation of N 2 (b l U u , b' 1 c\ l XT ) excited molecules in the band 98.0- 

102.5 nm. This source term was discussed and deduced in Ref. [26] and its expression for 21) 
case has been obtained in Refs. [19,20] 

T * -X x 

— I f 

x 0 

(2.14) 

p„ = ■•:,)' +(//-;'/„)' • 

Here Q h = 1.4x10 lf cm 2 is an average over the absorption band photo-ionization cross 
section of CE molecules, N u — 0 . 22 N is O 2 molecules concentration in air, a M 0.07 is a 
ratio of excitation rate of NS molecules by electron impact to their ionization 
rate, r * / r J{ = (25.33N / A r 0 -|- 1 ) 1 is a ratio of the total lifetime of excited molecule 

N 2 (b l U . ES c\ E*) to its radiative lifetime, and h s/ , ~ 0.1 mm is the characteristic size of 
streamer head. 

For dimensionless presentation of the results we use the sealing parameters: /7o 0.82x10 12 em 1 

for charged particle densities, £0 = 40.35 kV/em for electric field, and 1.091 kV for 
potential. The ionization threshold electric field in atmospheric air is £,/, 32.28 kV/cm. 

2 1 lb Initial and boundary conditions for transport equations 

Because the diffusion flux of ions is 4-5 orders of magnitude smaller than their drift flux, it can be 
neglected everywhere, even at the discharge front. lienee the electron transport equation is 
parabolic, and the ion transport equations become hyperbolic without diffusion flux. In accord to 
the type of equations, the boundary conditions read 

n — n t — n — 0 , if v -> x, or v — » ± x, (2.15) 

for charged particle concentrations at infinity; 

n n 0, if y = h , x < 0, (2.16) 

for positive and negative ion concentrations at the electrode surface; 

n = 0 , if y = 0 , v > 0 , and £ v > 0 , (2.17) 

for positive ion concentrations at the dielectric surface, when the ion flux is directed outward the 

surface. For negative ion concentration the similar condition reads 

n 0 , if y - 0 , x > 0 , and E y < 0 . (2.18) 


s,,h( x >y) = 


27TT,, 


) <;x P (- Q,.h S () ,p„) arc t a 11 


0 

•Cady,, 


P„ 
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For electrons, we use a general boundary condition coining from ab initio physics and valid for 
any surface. This condition expresses the hydrodynamic flux near the surface via the kinetic flux 
and reads 


r r 0(D ' U ) ( f 

I\ n t -= -alh 

’ • V v ' 

dt/ 


n V, 


— 4 - f t , for v = 0 , x > 0 and y — h t , x < 0 , 


(2.19) 


where 


o(I\) 


(1 


(1 - r)exp 


I< EX 
, T 


for E u < 0 


exp 

KEX 

+ 4 

1 - exp 

KEX)) 

, T ; 

T JJ 





L 


- 7 ,J , if y = h ,, x < 0, E 9 < 0 
0 if y = 0 , x > 0 


( 2 . 20 ) 


( 2 . 21 ) 


Here Vj - electron thermal velocity, A = electron mean free path, r =- surface reflection 
coefficient for electrons (hereafter it is assumed that r = 0 ), = elcetron secondary emission 

coefficient, _ electron flux from the surface {f oM ~ 0 for the dielectric surface and the anode 
electrode, and / = - 7 for the cathode electrode). 


Ihe electron secondary emission coefficient y s depends on electrode material and applied 
voltage. Usually it is of the order of 0.01 0.1 [27]. We assume that y s - 0 , because our 

computations showed that results do not notably depend on this parameter. 

The balance equations for the surface density of electrons, <r ( ,, and ions, a , read 

do do j 

Ot c)t 

The initial conditions are 


v = 11 - n . ii — 0 . o t (;r. 0 ) — (j\ 0 ) - 0 at / 0 , (2.23) 

where n in ~ background concentration of electrons and positive ions. Hereafter we use n m = 
10 cm , This parameter affects only the initial delay time of the discharge formation and does 
not affect the discharge evolution. 

2. / ,Ic Numerical algorithm and refinement of ionization source 

Because of a steep electric field and clectron-ion density gradients relevant to ionization wave 
front in the case of streamer formation and to cathode layer structure in the ease of negative 
exposed electrode polarity, the spatial size of computational cell should be less than 0.002 mm. 
This mesh was used for calculations discussed hereafter. Transport equations were solved using 
an implicit numerical algorithm with “upwind'’ approximation for the electron drift flux and 
Gauss elimination (sweep) technique for the electron diffusion. To resolve a strong nonlinearity 
due to the ionization source in the right-hand side of Fqs. ( 2 . 1 ) and ( 2 . 3 ), an iteration procedure 
was applied. Relative accuracy of iterations was 0.01. The Poisson equation for the self- 
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consistent electric field potential was solved using the improved Gauss-Seidel (upper relaxation) 
method with relative accuracy 10' 

If the £-field is strongly inhomogeneous in space, then the ionization rate constant should he 
obtained by solving the inhomogeneous Boltzmann kinetic equation. Nevertheless, a reasonable 
estimate of this rate constant can be made using the energy conservation equation for electrons. 
Our theoretical analysis gives the following expression for the ionization rate constant, which is 
corrected to account for the non-loeal ionization. 


k 


ul 


= k 


F:V{D t n,y 

Kn r E~ 


(2.24) 


The new model comprising the ionization rate constant (2.24) allow us to simulate the discharge 
evolution from the beginning to the time instant when the discharge stops propagating along the 
dielectric surface. This was demonstrated for both negative and positive electrode polarities. 


2. / Id Results and discussion 


The physical and numerical model for SDBD in atmospheric air has been developed accounting 
for air photoionization by UV radiation from discharge zone, approximation of non-local 
ionization source and universal boundary conditions on the electrode and dielectric surfaces. The 
refined model and new computational code allow for simulations of the total SDBD cycle 
including both the discharge formation and the discharge relaxation phase. 



t = 3 56 ns 
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Figure 3 Electron density ne/nO contours evolution at 
V = - 4.5 kV, rectangular region denotes the exposed 
electrode. 
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Figure 4 Streamer evolution predicted by the 
developed model; electron density nc /n0 contours 
for V +4.2 kV. 


The discharge formation phase lasts few tens of nanoseconds and creates plasma. The discharge 
evolves as a streamer for positive electrode polarity (Fig. 4) and as a diffusive discharge with a 
narrow cathode region for negative electrode polarity (Fig. 3). The predicted values of the 
discharge length and the surface charge density agree well with available experimental data for 
positive and negative electrode polarities (Figs. 5 and 6). This validates the developed physical 
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model and computational code. Although the discharge formation phase gives negligible 
contribution to the momentum and heat sources relevant to (low control applications, it provides 
the initial conditions for the relaxation phase. 

T- - -1 



Figure 5 Surface charge density a in 
units of nC/cm 2 for simulated streamer 
evolutional) -t 7.12 ns, (2)-t 14.2 
ns, (3) experiment [16]; V- +4.2 kV. 



Figure 6 Experimentally observed S1)B1) length [16]: ( 1 ) 
constant voltage U > 0 , (2) alternating voltage, (3) 
constant voltage U <i). Symbols show predictions of the 
developed model at V- +4.2 kV and l -5 kV (V = U ). 


The relaxation phase begins when plasma is redistributed in such a way that it shields the 
external electric field and reduces the air ionization rate in the discharge region to very small 
value. The ion-clcctron and ion-ion recombination along with the ion drift motion are main 
processes in the relaxation phase which lasts few microseconds. This phase effectively 
contributes to the momentum and heat sources relevant to flow control applications. 


02 , 

[ 0 
nnL ' 



°0Q, 

*.o$l 



(D 

0.0 r 

00 0.2 0,4 

0.6 0.8 

1,0 

12 

1.4 1.6 

1.8 

2.0 

02 L ^ 






(2) 



VOlo 



n n L ^ 

+1.040 

+>-050 




U.U 

O.O 0,2 0.4 

0.6 0.8 

10 

12 

^- 

1.4 1.6 

1.8 

2,0 


0,2 


o.o 


jl. 


+3.010 0.000 

~° 0 i 5 +J.G10 


(3) 


| 02 


0.0 0.2 0 4 0.6 0 8 1.0 1.2 1 4 1.6 1 8 20 

0,000 (4) 


0 000 - 0,001 

- -0.005 


0.0 




+?005 


Q0f 




VOiq 43.010 o0 o1 


0 0 0.2 0.4 0,6 0,8 1.0 1.2 1 4 1.6 1.6 2,0 


x, mm 


Figure 7 Contours of time averaged momentum source x-componcnt in 
units of 107 N/m3 for SDBD relaxation at V - -5kV: (1) t = 0.18, (2) 

0.62, (3)- 1.1, (4)- 1.8 ps. 

The spatial distribution of charged particles and their composition are quite different in the 
discharge formation phases and the discharge relaxation phase. During the discharge formation. 
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plasma occupies a thin near-surface layer of %0.05 mm thickness for positive electrode polarity 
(streamer discharge) and %0.02 mm thickness for negative electrode polarity (diffusive 
discharge). In the relaxation phase, plasma spreads outward the dielectric surface, and the 
thickness of momentum-source layer increases to %0.3 mm at both positive and negative 
electrode polarities. For a case of streamer relaxation at decreasing applied \oltage. the thickness 
of heat-source layer is approximately 2 times smaller (^0.1 mm), and for negative electrode 
polarity it is close to the thickness of momentum-source layer. 

The momentum-source distributions have a complicated structure in space. For positi\e electrode 
polarity, the length of positive force region is approximately 0.3-1 0 mm. For negati\c electrode 
polarity, the body force is an order of magnitude greater than in the case of positive applied 
voltage (its absolute value achieves 10 s N/m 3 as shown in Fig. 7). 

2.7.2 Accomplishments of the 2010 effort 

The model developed in 2007-2008 allows us to simulate the building up and relaxation of a 
single microdischarge only. The obtained data arc not sufficient to predict the momentum and 
heat sources induced by actual DBD-plasma actuators. These sources result from few tens of 
consequent microdischarges developed during each half-period of alternating applied voltage. 
Because numerical modeling of this bulk of microdischarges is extremely time-consuming, the 
problem has not been solved yet except for the ease of negative ramp voltage [6]. 

In 2010 we tried to estimate the momentum and heat sources induced by SDBI) using a 
phenomenological model of discharge action on the quiescent air based on existing experimental 
data and achieved understanding of SDBI) physics [28,29]. The primary objective of this work is 
to obtain analytical approximations of the body force and discharge length in order to recognize 
main trends of their dependences on applied voltage parameters and dielectric material properties 
and understand how to improve the SDBD flow control performance. 


2.1.2a Phenomenological model of SDBD induced body force 

The body force generated by a surface dielectric barrier discharge is due to momentum transfer 
from electrons and ions to neutral gas molecules [30] 

F — c(ii fl — 77,, — 7i t )E (2.25) 

Here e = electron charge, n r yn }) yn n ~ densities of electrons, positive and negative ions, 

respectively, and E - electric field strength. According to Hq. (2.25) the body force is 
proportional to sclf-consistent E-field value and positive or negative charge density surplus. 

Electrons and ions are produced during the building up of microdischargc, which lasts few tens 
of nanoseconds. Then, between the consequent microdischarges, they recombine, convert to 
another type of ions and/or move to the dielectric surface leaving the gas volume. This relaxation 
phase is characterized by a microsecond time scale. The electrical charge generated in a 
mierodiseharge due to air ionization is divided into three parts. The charges of sign that is 
opposite to cleetrode potential sign move to electrode and form conductivity current in the 
external electrical circuit. The charges of the same sign as the sign of electrode potential drift out 
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of the electrode and form both the surface charge layer on the dieleetnc surface and volumetric 
eharge. These charges screen the external eleetrie field and lead to the microdiseharge end off. 

Inside a microdischarge series, which is relevant to the voltage half-eyele of constant polarity, 
the eharges having the opposite to electrode polarity sign are not delivered to the surfaec. 
Accordingly, no charge recombination oecurs on the surfaee during the microdischarge series, 
and the surface-charge relaxation time is mueh greater than the half period of applied voltage. 
This means that the recovery of breakdown condition and repetition of mierodiseharges in a 
mieroseeond seale observed in experiments is possible only due to the increase of applied 
voltage and not due to the decrease of surface charge. 


Numerical simulations [28] of a set of mierodiseharges at a ramp voltage showed that at the 
beginning of the mierodiseharge series the time interval between the neighbor mierodiseharges is 


inversely proportional to the voltage slope; Ar mj , - I —) ~ —, where I 0 = voltage amplitude 


dV 


dt 


1 




and f v = voltage frequency. The number of mierodiseharges in a series is 

TAK-K,) 


N 


2 A Tjr a 




(2.26) 


where V h<{ = breakdown voltage for the first mierodiseharge, and T — / 1 is period of applied 
voltage. 


The number of mierodiseharges is proportional to the difference between the voltage amplitude 
and the breakdown voltage and may be determined as 


A 


„(*) = 


m - v;„ 

sv 


(2.27) 


where 6V is the voltage increase, whieh is neeessary to reeover the breakdown condition for 

T <M 

igniting a new mierodiseharge. The relation between 6\ and At tl is At, —-. 


The volumetric charge created in a mierodiseharge has an opposite to eleetrode potential sign 
Therefore, the body force, which is averaged over time and integrated over spaee, is positive 
(directed out of the electrode edge) both for positive and negative eleetrode polarity. This was 
shown by straight and time resolved body force measurements [31-34]. Following the eommonly 
used notation [31] we call the discharge at a negative going potential half cycle as a forward 
stroke and the discharge at a positive going potential half cycle as a backward stroke. 

The volumetric eharge surplus leaves the volume due to drift onto the dieleetrie surfaee in 
tangential direction (along the dieleetrie surface), because the normal component of eleetrie field 
has been screened by surfaee eharge at the stage of mierodiseharge building up. The volumetrie 
eharge in air is primarily presented by molecular ions. In the ease of forward stroke, when 
electrons start to form the volumetrie eharge, it is due to rapid electron attachment to oxygen 
molecules and rapid conversion of atomic ions 0 into molceular ions 0 5 [20]. In the ease of 
backward stroke, the molecular ion structure of volumetrie eharge is natural because positive 
molceular ions arc created due to air ionization and further they just ehange their composition in 
ion conversion reactions [20]. Hence the relaxation time of the volumetric eharge may be 
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estimated as the drift time of molecular ions along the dieleetrie surfaee: A t j « L / v %ir , where 
= ion drift veloeity. According to the numerical simulation [29] the average E-field for 

mierodiseharge decaying air plasma is around 5 kV/em, the ion mobility approximately equals 2- 
3 cm7(Vs). This gives % 10 em/s. 

The discharge length L depends on the applied voltage amplitude Jo* the dielectric thickness ci 
and the dielectric relative permittivity e [16]. As a reference point we consider l u « lOkV, 
(I % 1 mm, e %5-10. For this conditions, L % 10 mm [16] and At « 100//.S. The experimentally 

observed time interval between microdiseharges is primarily governed by the voltage increase, 
and it is Ar i(i</ % 15 ps. Accordingly, the volumetric charge relaxation time At is much greater 

than the time interval between microdiseharges, Ar A t wiI . Relaxation of the volumetric 

charge mainly occurs after the end of a microdischarge series, when the voltage slope changes its 
sign, Henee the volumetric eharge aeeumulates during the development of microdischarge series, 
and the body foree indueed by this increasing eharge grows during the voltage half cycle. For the 
forward stroke, accumulation of the negative volumetric charge has been modeled numerically in 
Ref. [28]. The force increase during the voltage half eyele has been demonstrated experimentally 
in Refs. [33,34] for forward and backward strokes. 

The accumulated volumetric charge is proportional to the number N,„j of microdiseharges in a 
series. For normal SDBD developments, this number remains constant for the constant voltage 
amplitude and does not depend on the voltage frequency (see Hq. (2.26)). In this case, a set of 
mierodiseharges ends off when the voltage achieves its amplitude value F 0 , and the voltage slope 
changes its sign [16,31]. 

The and corresponding surfaee and volumetric charges increase with J o until they reach 
saturation levels. We suppose that this occurs when the volumetric charge, which is accumulated 
in a mierodiseharge series, beeomes high enough to contribute into the external voltage screening 
in the same extend as the surface charge. If the saturation is achieved, the set of microdiseharges 
ends off before the voltage reaches its amplitude value. The corresponding current signal of 
mierodiseharge series was observed experimentally [35]. The maximal body force was observed 
for the maximal volumetric charge, and N tm { was reaehed at the saturation point 

Because of great difference between the characteristic gasdynamic time 
T q ~ L / (7 X ~ 10 ’* — 10 l s and the mierodiseharge evolution time - 10‘ {1 s (L is the discharge 
length, and (7 X is the frcc-strcam velocity), the momentum and heat sources averaged over the 
alternating voltage eyele are relevant to aerodynamic forcing. The body foree, averaged over the 
voltage cycle, includes integration of the instantaneous foree over the forward and backward 
strokes 

/ 2 7; 

/" = — f F ( t)dt +— f F (l)ilt . 

7 J j J 

According to theoretical predictions for a single mierodiseharge evolution [29], the backward 
stroke foree F (t) should be an order of magnitude smaller than the forward stroke force/’(/). 

This is because of quite different plasma and clcctrie field spatial distributions generated at 
positive and negative electrode polarity (see Section 2.1 I). In the experiments [33,34], the body 
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force were approximately the same for the forward and backward strokes. These measurements 
provided an average (over the exposed electrode length) force. The recent experiment [32] 
showed that the body-force spatial distribution along the electrode length is extremely 
inhomogeneous. Generally, the body force for the negative going voltage half-cycle is much 
greater than that for the positive going voltage half-cycle. But at some points on the electrode 
edge the reverse ratio takes place, boregoing means that the correct estimation of the body force 
value could be made by accounting for the forward stroke force only. The backward stroke force 
contribution may either be negligible or increase the result not more than by factor 2. 


2.1.2 b Body force estimation Jor negative going voltage 


Consider the SDBD evolution at a normal regime far from the saturation limit. For negative 
electrode polarity, the microdischargc process is similar to that observed in a glow discharge 
[29]. The cathode layer is formed after air breakdown in the vicinity of the exposed electrode 
edge. Electron-ion generation due to air ionization occurs primarily inside this cathode layer. The 
ions are moved to the cathode by strong electric field, and electrons drift in the opposite direction 
towards the dielectric surface. Assume that each microdischarge generates approximately the 
same amount of electrons possessing a charge A<y. The main part of these electrons rapidly 
leaves the discharge volume reaching the dielectric surface, and the rest part is converted into 
negative 0 2 ions due to attachment to O, molecules [20] creating a volumetric charge Ar/ of 
negative ions. 


The probability of the negative ion formation on electron drift way to dielectric surface is 

»J(t) 


estimated as pft) 


, where is til ~ attachment frequency of electron to 0 molecule, /(/.)- 


KEAt) 

current value of the surface charge length (discharge length), A electron mobility, and E ,(/)" 
current value of electric field tangential component above the surface charge layer. Electrons 
drift in the tangential direction because the nonnal electric field component equals to zero due to 
the surface charge screening. The surface charge layer length is 


XmlAqE ft) 

eVM) 


L — nia x(/(/)), 


(2.28) 


where V= average potential of the dielectric surface at a surface charge location. The negative 

HO 


ion residence time inside a volume is estimated as At. 


E Eft) 

mobility. 

Estimation of the time averaged body force inside a microdischarge series reads 


, where I\ ~ negative ion 


- 1 V 

h % T ~J ( E( f )E jt ( f )dE 


(2.29) 


Following available experimental data we assume that A r iW (/) and /(/) are proportional to /, the 
current value of accumulated volumetric charge is expressed as 


16 






( 2 . 30 ) 


Then 


q tt (i) *-N mil {t) -Aq - p„(t). 




Integration in Eq.(2.31) finally gives the body force 

3KK.E. 


(2.31) 


(2.32) 


F-quation (2.32) is valid for / Ar w <1/2. More correct expression, which accounts for the body force 

decrease for / Ar n > 1 /2, is 


- ^ vJ.S iq exp ( _ /ATn )[l _ exp (-/At )j. (2.33) 

oft 

The experiments [32] show that all mierodiseharges in a series are started from the exposed 
electrode edge and not from the edge of the surface charge layer created by previous 
mierodiseharges . Hence, in order to estimate A q,N wti and /., we have considered an asymptotic 

analytical solution for eleetrie field near the electrode edge. The air breakdown condition reads 
[27] 

j*a T (r)dr = ln(l + 7 '), (2.34) 

c, 

where the Townsend coefficient is integrated along the line of £-field Ce from the electrode 
edge to the dielectric surface, and 7 = coefficient of secondary electron emission from the 
cathode surface. A semi-empirical expression of a T for the air is a T /lpexp( Bp/E), 
where A- 15 em ’Torr \ B - 365 Vcm 1 forr , and p ~ air pressure in Torr [27]. 

The breakdown condition (2.34), with substitution of the aforementioned E-field near the 
electrode edge after appropriate analysis and calculations, gives the breakdown voltage 


7T B 


pd ln(l f 7 ') 


2A 


(2.35) 


The E-field contours Ce satisfying the condition (2.34) are ended on the dielectric surface at 
some distance from the electrode edge. Hence the surface charge layer, which is created by 
drifting to dielectric surface electrons, does not flank to the electrode edge. Numerical simulation 
shows that the distance l hfl between the electrode edge and the starting point of surface-charge 

layer is about dll Hereafter we use 


To get the voltage increase^ 7 , which is necessary to compensate the screening effect of the 
surface charge and restore the breakdown condition, we have analytically estimated the electric 
field inside the gap between the electrode edge and the surface charge layer. By variation of the 
air breakdown condition in the gap we obtain 
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5V(kV) « 2d(crn), kV. (2.36) 

Using the foregoing estimates of t>V % V M , E x i , N mti and Ai/, the discharge length and the body 
force are evaluated as 


L{ cm) « 0.1 • V;,(kV) 1-V • 

K) J 

(2.37) 

F(K/m )« 10 7/ (kH ^lV(kV)(l 1M 1 . 

d{cm) ( , 

(2.38) 


Figure 8 shows that the theoretical estimate (2.38) agrees satisfactorily with the experimental 
data [36] for dielectrie layers of different materials and thiekness. Some disagreement for near 
breakdown voltages can be attributed to the inaccuracy of the breakdown voltage estimation 
because of inaeeuracy of the used asymptotie /T-field value near the eleetrode edge. 



Figure 8 Theoretical estimations of the body force (lines) against the experimental data [36] (symbols) 
for 1,1a - teflon (e= 2), cl 6.35mm,/.=2.1 kHz; 2,2a teflon, cl 3.18 mm, j\ 2 kHz: 3,3a kapton 
(c =3.9), d =0.15 mm,/ 4.4kHz. 


2.1.2c Discussion 

The developed phenomenological model sheds light on physical processes governing the ignition 
and ceasing of mierodischarges for the negative going half-cycle of applied alternating voltage. 
This model clarifies reasons of the body force saturation observed in experiments. The obtained 
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analytical formulae capture basic dependencies of the body force and discharge length on the 
applied voltage parameters, dielectric thickness d and its relative permittivity e . 'I he model is 
validated against the experimental data providing the body force variation versus voltage 
amplitude for different dielectric materials and their thicknesses. 

For negative electrode polarity, a gap between the exposed electrode edge and the surface charge 
layer provides an easy start for consequent microdischarges. Because the electric field inside this 
gap remains rather high, even when a mierodiseharge ends off and relaxes, a small increase of 
applied voltage is sufficient to restore the air breakdown condition and start a new 
mierodiseharge. This voltage addition primarily depends on the dielectric thickness </, while it is 
weakly affected by the voltage value and 5 . The number of mierodiseharges during the voltage 
half-cycle is proportional to the voltage amplitude and does not depend on the voltage frequency. 
Hence the time interval between mierodiseharges is inversely proportional to the voltage 
amplitude and its frequency. 

For negative going voltage, the discharge length is proportional to the voltage amplitude and 
practically does not depend on d and 5 . This result looks surprising because the experimentally 
observed discharge length for alternating voltage decreases with £ [16] This could be explained 
by the fact that the discharge length is much greater for the positive going voltage half-cycle. 
Accordingly, for alternating voltage, its dependence on £ is associated with the discharge length 
at positive going voltage. Additional experiments are needed to clarify this issue. 

The proposed model captures experimental dependencies of the body force on voltage at 
different d and £. This indicates that the volumetric force is primarily induced by the negative 
going voltage half-cycle. The force associated with the positive going voltage half-cycle should 
not change the result significantly. 

The body force is due to accumulation of the volumetric negative charge carried by negative 
long-lived 0 2 and 0 :i ions. This accumulation is proportional to L y giving a strong dependence 
of the body force on applied voltage. The probability of negative ion formation due to electron 
attachment is proportional to l: and the residence time ol negative ion in the volume is also 
proportional to L~. 

To increase the force and improve the SDBD actuator performance, one has to increase the 
volumetric charge of negative ions. This could be achieved by artificial increasing of the 
discharge length. The latter can be done by shifting the encapsulated electrode downstream of 
the exposed electrode or by using a dielectric of small conductivity. 

The idea to use a “conducting” dielectric was proposed and tested in Princeton two years ago 
[37]. They could not improve the actuator performance. The reason of this failure may be due to 
the fact that they have used nanosecond pulses instead of a sinusoidal voltage of moderate 
frequency. They tried to decrease the surface charge in the time scale between the pulses 
preventing the negative ion accumulation. To get a positive effect one should use a sinusoidal 
voltage (supporting a mierodiseharge series and negative ion accumulation) and supply the time 
of surface charge reduction due to finite dielectric conductivity close to the voltage period. For 
the material used in experiment [37] (linen based phenolic) this time corresponds to frequency of 
applied voltage -1 kllz . It is expected that further studies of this problem can lead to significant 
improvements of the SDBD performance. 
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2.2 SDBD-induccd momentum and heat sources 


Using the body-force absolute values predicted by the developed phenomenological model and 
numerical results for the spatial distributions of momentum and heat sources, we obtain 
analytical approximations of the volumetric force components (F r , F ) and the heat source Q 


F T {x.y.l) = if(t)F t ,,{x,y), 

F„(x.y-1) = ■p(t)F ild (x,y)> 

Q(x,y,l) = <p(t)QAx s t), 

where 

K, - AM-r)f(y)< F„ = Q, - ■ 


The functions g tl (x ), f(y) and f (y) are expressed as 


a A*) 



. < X < X ,, 


/('/) = AjV exp 

/,(?/) = A,?/exp 




(// — Ai)’ 


,A 


= 7547 m \ 
= 3-10' m 


a f = 10 "?// 


'.a = 4 • \if m 


(2.39) 


where l tl — 5 mm is the discharge length, //, = 0.1mm and — 0.005 mm are vertical sizes 
of the discharge region, .r, = 2.5 em, x. 2 = x , T /,. The ^(/) is a step-function ^(1) = 0 for 
I < 0, ^»(/) = 1 for / > 0. The constants in Hq. (2.39) arc 

A j , l = ±Ur N/m\A wl = -2-U)' N/nA.A* = 10 7 M/m’. 

Plus (or minus) in the expression for/l jr/ correspond to the upstream (or downstream) shift of the 
exposed eleetrode with respect to the buried electrode. 


2.3 Boundary-layer flow control using SDBD actuator 

To shed light on aerodynamic effeets produced by SDBD actuator, wc consider a relatively 
simple flow - the laminar boundary layer on a Hat plate in a free stream of low 7 subsonic speed. 
The plate has zero thickness and sharp leading and trailing edges. The computational domain and 
Cartesian coordinate system arc shown in Fig. 9. Part of the bottom boundary corresponds to the 
plate surface 0 < x < 0.1 in. The inlet boundary is located at x = 0.05 m. and the outlet 

boundary is located at i = 0.2rn. The 2D SDBD forcing is produced in the region 

0.025 in < x < 0.031 m. The computational grid has approximately 210,000 hcxahedral cells: 
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300 grid nodes are in the boundary-layer region in the vertical direction, and 580 grid nodes are 
along the plate surface. Clustering of the grid nodes is used in the \ieinity of SDBD actuator, 
0.020 m <x< 0.050 m. 


1- 

Top 



Inlet 



Outlet 

V ‘ 

SDBD Region 

1 1 , 



Symmetry 

Plate LE v 

Plate TE j Symmetry 



Figure 9 Computational domain and coordinate system 

Although the flow has small subsonic speeds, the gas compressibility is taken into aecount to 
simulate effects of the SDBD-indueed heat source. 2D unsteady Navier-Stokes equations are 
integrated using an implicit second-order finite-volume method. The governing equations are 
approximated by a conservative scheme. The viscosity-temperature dependence is calculated 
using the Sutherland law with constant 1 10.4K. The fluid is a perfeet gas having the specific heat 
ratio 7 = 1.4 and Prandtl number Pr = 0.72. 

The freestream parameters are: velocity i/ x = 10 m/s, pressure 10 Pa; density 

p x =1.2 kg/m 3 , temperature 290 K, the Reynolds number is 

Re x — — 1.73 x K) 1 , where x x - 0.025 m is the SDBD leading edge coordinate, 

which is chosen so that the boundary layer thickness is <*>(.r,) ^ 1 mm. 

The freestream conditions are imposed on the inlet and top boundaries. On the outlet boundary, a 
linear extrapolation of dependent variables is used. The boundary conditions on the bottom 
boundary are: the flow symmetry condition upstream from the plate leading edge and 
downstream from the plate trailing edge; no-slip condition and the adiabatie wall temperature on 
the plate surface. 

2.3.1 Steady SDBD forcing 

In this series of computations, the source terms are initiated at the time instance f 0 and kept 
constant for / > 0 . The unsteady aerodynamic problem is solved with the time step 5* 10 u s. 

Consider the case when the SDBD momentum source is directed downstream. I he transversal 
component of vortieity field (sealed from 3 10 1 s l to 3-10 1 s' 1 ) and temperature field 

(sealed from 290 K to 291 K) are shown in Fig. 10 at different time instants. The SDBD 
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actuator generates a longitudinal near-wall jet. The jet head floats to the upper boundary-layer 
edge and rolls up to a vortex propagating downstream over the boundary layer. 

/ - 3-10 1 s 



Figure 10 Vorticity and temperature fields at different time instants, steady downstream SDBD forcing. 

Another ease corresponds to SDBD acting in the upstream direction. Snapshots of the transversal 
component of vorticity Field (sealed from 3 • 10 1 s to 3 -10 1 s 1 ) and the temperature field 

(sealed from 290 I\ to 291 K ) are shown in Fig 11. T he SDBD actuator pushes the near-wall 

fluid upstream and generates a counter-flow jet. This jet interacts with the boundary-layer flow 
that leads to excitation of strong concentrated vortices propagating downstream. In this ease, the 
SDBD actuator works as a vortex generator. 
t = 3-10 1 s 






















/ = 7.5 • 10 ' s 




l = 15 10 1 a 




t. = 50-10 1 s 




Figure 11 Vortieity and temperature at different time instants, steady upstream SDBD forcing. 

2.3.2 Unsteady SDBD forcing 

An unsteady SDBD forcing is simulated by repetition of SDBD cycles as schematically shown in 
Fig. 12. It is assumed that when the discharge is on, the momentum and heat sourees are 
described by the formulae of Section 2.2. When the discharge is off, the sourees are zero. The 
momentum and heat sourees are modulated with tune as 


v( 0 = 


sin u?/,sin jjl > 0 
0 , sin u)i. < 0 


(2.40) 


Here the circular frequency is j - 2ttU x / L f ss 10 1 s 1 that corresponds to the modulation 
period T = 6 * 10 1 s. 


i 
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F igure 12 Schematic distribution of SDBD induced sourec modulation. 
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Snapshots of the transversal vortieity field (scaled from -3*10 s to 3*10 1 s' 1 ) and the 

temperature field (scaled from 290 K to 291 K ) are shown in Fig. 13 for the case of 

downstream forcing. The SDBD actuator generates a sequence of longitudinal near-wall 
pulsating jets. In the far field, these jets roll up to vortices propagating downstream over the 
upper boundary-layer edge. 



Figure 13 Vortieity and temperature fields at different time instants, unsteady SDBD acts 
downstream. 

The case of upstream SDBD forcing is illustrated by snapshots of the transversal vortieity field 
(scaled from -3*10 1 s' 1 to 3*10 1 s 1 ) and the temperature field (scaled from 290 K to 

291 I\ ), which are shown in Fig. 14. The SDBD actuator generates strong concentrated vortices 

which propagate over the boundary layer and strongly perturb the near-wall How. 

The foregoing numerical solutions demonstrate that by choosing appropriate distributions of 
applied voltage it is feasible to achieve quite different aerodynamic effects on the near-wall flow. 
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This flexibility and very short characteristic timescales make SDBD actuators attractive for 
control of the near-wall flows associated with the boundary-layer separation, laminar-turbulent 
transition or a combination of both. 

L = 5 -10 1 s 



t = 10 ■ 10 “' 5 



f = 50-10 ' s 



Figure 14 Vorticity and temperature fields at different time instants, unsteady SDBD acts upstream. 


2.4 Vortex flow over a supersonic delta wing at angles of attack 

Data relevant to vortex structures and their control for delta wings in supersonic flow are very 
limited. This motivated us to perform numerical simulations of the vortex flow past a delta wing 
with sharp leading edges at supersonic speeds (freestream Mach number 1.5) and estimate 
feasibility of the vortex flow control using SDBD actuators. 

2.4.1 Problem formulation and numerical method 

Consider 3D viscous flow past a delta wing at the free-stream Mach number 1.5 and various 
angles of attack. The delta wing configuration is shown in Figure 15: the apex angle is 60°, and 
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the leading-edge sweep angle A — G0°. The wing has zero thickness and sharp leading and 
trailing edges. 

Numerical solutions are obtained using an implicit finite-volume method. Three-dimensional 
Navier-Stokes equations are approximated by a conservative scheme. The flux vector is 
evaluated by an upwind, flux-difference splitting of Roe [39]. MUSCL algorithm is applied with 
the third order TVD space discretization [40]. An Euler implicit discretization in time of the 
governing equations is combined with a Newton-type linearization of the fluxes to obtain the 
system of algebraic equations [41]. This system is solved using a point Gauss-Seidel scheme. 

The no-slip boundary conditions are imposed on the delta wing surface. Temperature on the 
wing surface equals to the adiabatic wall temperature. On the outflow boundary, the unknown 
variables are extrapolated using the linear approximation. On the inflow boundaries, the 
conditions correspond to free stream. Due to flow symmetry, problem is solved only for one half 
of the delta wing with symmetry conditions on the plane y — 0. 



Figure 15 Delta wing configuration and coordinate system. 


2.4.2 Cross-flow topology 

In order to specify the parameter space for computations, wc have analyzed available data on 
vortex flows past supersonic delta wings. Different cross-flow patterns, which were observed 
experimentally and numerically over the leeside of delta wings at various angles of attack [42- 
44], arc summarized in Fig. 16 in terms of Mach number normal to the leading edge 

M s = M x «>sAJ(l + sin 5 a taif A) and the angle of attack normal to the leading edge 

n N = tan 1 (taiitt/eosA). With the help of the flow visualization data (vapor-screen 

photographs, tuft and oil-flow photographs) and surface pressure data, the flow patterns were 
subdivided into seven types: 1) classical vortex, 2) vortex with shock, 3) separation bubble with 
no shock, 4) separation bubble with shock, 5) shock with no separation, 6) shock induced 
separation, 7) no shock and no separation. 

More detailed topology suggested in Ref. [44] is shown in Fig. 17. The gray lines (1, 11 etc.) 
depict boundaries shown in Fig. 16. The colored lines 1, 2, 3, 4 and 5 show the new details: the 
green line 1 separates the flow regimes without and with the secondary separation; the yellow' 
line 2 separates the flow regimes with tertiary separation; the region above the blue line 3 
corresponds to the flow with a shock wave under the primary vortex; the region below the blue 
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line 3 corresponds to the flow with no shock under the primary vortex; the brown region 4 
corresponds to transition between regimes with and without shock wave above the primary 
vortex; the red line 5 separates the flow' regimes with and without shock between the primary 
vortices. 


O Classical vortex 
CJ Separation bubble 
O S^ock-induced bubblp 



Figure 16 Cross-flow topology reported in Ref. [42] 
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Figure 17 Refined cross-flow topology suggested in Ref. [44], red line with squares our points. 
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In all eases, the flow pattern is symmetric heeause of sharp leading edges of a thin delta wing. 
Note that the cross-flow diagrams are presented in terms of inviscid parameters o v and M s , 
while flow separations depend on the Reynolds number. Therefore, the flow regimes with 
secondary and tertiary separation may disagree with these diagrams. The number of separations 
ean increase as the flow approaches the trailing edge [43]. 

Using the aforementioned data we specify the following set of flow parameters for our 
computations: the freestream Maeh number M x — 1.5, the stagnation temperature 

T* = 300 K, the Reynolds number Re x — p\U\L'j — 2 x l() f \ where U 1 m is length of 

the delta wing; and angles of attaek: o = 0°,5°, 10°, 15°,20°,25°,30°. The corresponding points 
in the M s — a v plane are shown in Fig. 17 by the red line with squares. 


2.4.2 Vortex flow patterns and aerodynamie forees 


Our numerical solutions gave the cross-flow patterns consistent with the diagrams shown in Fig 
17. As an example. Figures 18-22 illustrate the case of a = 30° that corresponds to the regime 
with shocks above and under the primary vortex. These shoeks are clearly seen in Fig. 18. The 
cross-flow streamline patterns (Fig. 19) show that the secondary separation is formed in the wing 
mid-span under the primary vortex. Substantial differences between the pressure distributions at 
the stations x = 0.25 and x — 0.5 (Fig. 20) indicate that the conical pattern breaks down 
somewhere in between of these stations. This is confirmed by 3D distributions of vorticity (Fig. 
21) and streamlines (Fig. 22). 



Figure 18 Total pressure ratio field at x — 0.25 (left) and x 0.5 (right), o — 30°. 
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Figure 20 Pressure coefficient for x = 0.25 (green line) and x = 0.5 (blue line), a - 30°. 



Figure 21 Isosurface of constant vorticity, a ~ 30°. 



Fig. 22 Streamlines near the delta wing, n = 30°. 
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The integral lift and drag coefficients are shown in Fig. 23 as functions of the angle of attack. 
The wing aerodynamic performance degrades as the angle of attack exceeds 25°. This correlates 
with the vortex breakdown that occurs near the wing trailing edge at a ~ 25° and quickly 
moves upstream as a increases. Detailed analysis of the flow patterns showed that the wing 
aerodynamic performance depends on the locus of the conical flow breakdown over the wing 
surface. 



2.5 Vortex flow control strategy 


As shown in Section 2.4, the delta wing aerodynamic performance depends on the locus where 
the conical vortex structure breaks down. This suggests that the high-angle-of-attack flow control 
strategy should be focused on the control of vortex breakdown over the wing lceside. This 
phenomenon is associated with a large-scale instability of vortices that is characterized by rapid 
deceleration and eventual stagnation of the vortex core flow as well as growth of the vortex core 
cross-section. Much effort has been attempted to understand the physics of this flow 
phenomenon (see the review articles by Hall [45], Leibovich [46], Escudier [47], Delery [48], 
Althaus et al. [49], Lucca-Negro and O'Doherty [50]). 

The need to control the vortex breakdown over a delta wing at high angles of attack has 
motivated researchers to propose various flow control techniques [51] including: mechanical 
devices such as leading-edge flaps [52,53], apex fences [54], canards, strakes, leading-edge- 
extensions (LFXs), or double-delta w ing [55-62] and pneumatic techniques such as trailing-edge 
blowing [63-68] or suction [69], along-the-core blowing technique [70-73], spanwise blowing 
[74,75] and some other blow mg/suction techniques [76,77] with different blowing/suction 
locations and orientations. 

It is generally believed that the trailing-edge blowing or suction techniques delay vortex 
breakdown by decreasing the downstream pressure gradient, and along-thc-core blowing 
technique delays vortex breakdown by increasing the axial velocity, i.e., a reduction in the swirl 
number. Dixon [74] believed that spanwise blowing on the wing provides sweep-like effects as 
the blowing jets are entrained in the leading-edge vortices. Bradley and Wray [76] credited the 
success of their blow ing technique to the increase in the vortex stability, which to some extent is 
related to the longitudinal flow in the vortex core. For canards, strakes, leading-edge-extensions, 
or double-delta wing [55-62], these devices manage vortex breakdown by modifying the flow' 
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field and inducing a nonunifonn distribution of local angles of attack at the wing, leading to the 
generation of a non-conical vortex formation and the corresponding delay in vortex breakdown. 

SDBD actuators may be advantageous compared to the aforementioned techniques. This 
motivated Visbal and Gaitondc [38] to conduct numerical simulation of the DBD actuator effect 
on the vortical flow above a swept delta wing at high angles of attack in low-speed (Mach~0.1) 
free stream. With the actuator located near the apex, significant movement of the vortex 
breakdown location and a dramatic transformation of the shear-layer sub-structurcs were 
predicted. Hereafter we discuss results of our numerical simulations of the SDBD effect on the 
vortex flow past a supersonic delta wing introduced in Section 2.4. 

2.6 Vortex flow control using symmetric SDBD forcing 

Consider the wing-apex SDBD, leading-edge SDBD and multi-element SDBD actuators 
schematically shown in Fig. 24. These configurations arc symmetric with respect to the wing 
centerline, and computations have been carried out for one half of the total flow region w ith the 
symmetry conditions imposed on the plane z — 0. 



t* 

(b) 


(a) 


Figure 24 Red lines shows SDBD regions for (a) the wing-apex actuator, (b) the leading-edge actuator 
and (c) the multi-element actuator. 

For the wing-apex and multi-element actuators, the body force components are 



(2.41) 


(2.42) 


For the leading-edge actuator, the body force F t is perpendicular to the leading edge. In this 
case, the streamwise and spanwise components are determined as / F { cos A, 
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F = F l sin A. The heat source tenn in the energy equation is approximated as Q — (y). 

For all cases, F 0 = 10 1 Njm *, y {] — 3x10 ’ ro, Q {] = 2xl0 9 Wjrn m \ the width of SDBD 
region is 1.5 cm. 


2.6.1 Solutions for the wing-apex actuator 


Consider the wing-apex actuator (Fig. 24a) with the boundaries of the SDBD forcing region 
r, = 3 cm and x 2 = 4.5 cm. By changing sign of F } wc have simulated SDBD acting 
downstream ( F L > 0) and upstream ( F t < 0). In the cases of a ~ 0° and a — 5°, there is no 
appreciable effect of SDBD on the flow field because there is no global separation from the wing 
leading edges. For o = 10° (Fig. 25), the SDBD actuator affects the wall temperature field and 
the streamline pattern. First evidence of the vortex breakdown is observed near the wing trailing 
edge in the case of no SDBD forcing (a). 



(a) SDBD is off 



For o — 20° (SDBD is off) a well-developed vortex is seen in a mid-chord station (Fig. 26a). 
The burst locus moves downstream with the SDBD acting downstream (Fig 26b), and it moves 
upstream with the SDBD acting upstream (Fig. 26c). This example demonstrates that the vortex 
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breakdown locus can be controlled by the SDBD forcing near the delta-wing apex. T his is 
consistent with the low-speed modeling of Visbal and Gaitonde [38]. However the integral 
aerodynamic forces (lift and drag coefficients shown in Figs. 27 and 28) arc weakly affected. 
Presumably the aerodynamic loads at supersonic speeds arc not so sensitive to the vortex burst 
locus. For a > 25°, the influence of DBD on the vortex breakdow n is not so clear, because the 
breakdown point is very close to the wdng apex in all three cases. 



(a) SDBD is off 



For sufficiently high angles of attack, the How field may evolve w ith time. To clarify this issue 
we performed direct numerical simulation of unsteady vortex fields for the ease of a = 20°. As 
expected the primary-vortex breakdown slowly moves along the wing surface. The secondary 
vortex reveals more unsteady oscillatory behavior. This unsteadiness is observed both with and 
without SDBD forcing. Nevertheless, appreciable migrations of the burst locus weakly affect the 
lift coefficient C { (it varies in the range 0.836<f r / <0.840). 
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Figure 27 Lift coefficient versus angle of attack. 



Figure 28 Drag coefficient versus angle of attack. 


34 































2.6.2 Solutions for the leading-edge actuator 


Consider the leading-edge SBD actuator shown in fig. 24b. The boundaries ofSDBD forcing are 
specified as a;, = 3 cm, x 2 — 6 cm and x. A — 90 cm. Figure 29 shows the drag polar, data for 
the wing-apex actuator are also presented for comparison. In this figure: the red line corresponds 
to the case of SDBD off; the blue line corresponds to the case when wing-apex SDBD induces 
the momentum source in the positive x-direction (downstream forcing); the green line 
corresponds to the case when the wing-apex SDBD induces the momentum source in the 
negative x-direction (upstream forcing); the magenta square corresponds to the case when the 
leading-edge SDBD induces the momentum perpendicular to the leading edge at a —20°. 
Additional computations were carried out to distinguish the SDBD momentum effect from the 
SDBD heating effect. The black triangle corresponds to the case when the momentum source is 
included into Navicr-Stokcs equations while the heat source is not The cyan triangle 
corresponds to the opposite situation the heat source is on while the momentum souree is off 
In all the cases, SDBD forcing produces small effeet on the wing aerodynamic performance. 

Figure 30 illustrates the influence of SDBD on 3D streamlines and the wall temperature at the 
angle of attack o = 20°. When SDBD is off (Fig. 26a), the vortex burst is observed in the mid¬ 
chord station. The momentum souree (Fig. 30a) and the heat souree (Fig. 30b) lead to an 
appreciable downstream shift of the vortex burst locus. Nevertheless, they weakly affect the 
integral aerodynamic forces shown in Fig. 29. 
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(a) (b) 

Figure 30 The leading-edge SDBD effect on the vortex flow field at o = 20 ; (a) lhe momentum source 
is on while the heat source is off; (h) the heat source is on while the momentum source is off. 


2.6.3 Solutions for the multi-clement actuator 

Consider the multi-element configuration comprising five SDBD regions (Fig. 24c). This 
configuration resembles the SDBD actuator used for the low-speed wind tunnel experiments in 
1TAM [78]. The coordinates of the SDBD-forcing boundaries, x iit and .r , are given in Table 2. 

Figure 31 shows the five-element SDBD effect on the flow streamlines and the wall temperature 
pattern. Similar to the configurations considered in previous sections, this actuator causes 
appreciable downstream shift of the vortex burst. Nevertheless, the integral aerodynamic 
coefficients vary in the range of 3% only. 


Table 2 


Number of SDBD strip 



1 

0.03 

0.045 

2 

0.15 

0.165 

3 

0.4 

0.415 

4 

0.6 

0.615 

5 

0.8 

0.815 
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(a) SDBD is off (b)SDBDison 

Figure 31 The five-element SDBD effect on streamlines and surface temperature, n - 20° . 


2.7 Vortex flow control using asymmetric SDBD forcing 

Herein we discuss numerical simulations of aerodynamic moments induced by asymmetric 
SDBD forcing on the aforementioned delta wing. The problem is formulated for the total delta 
wing configuration without applying the symmetry conditions on the symmetry plane. This 
allows us to investigate feasibility of the rolling moment control by asymmetric SDBD forcing. 

The computational grid is obtained by doubling of the grid used for the symmetrical problem. It 
has approximately 9 x 10 6 nodes. In the boundary-layer and leading-edge regions, the grid nodes 
arc clustered to resolve fine flow structures. 



Figure 32 Boundary layer profiles. 
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First series of computations was conducted for the case of n — 5°. Figure 32 shows the 
boundary layer profiles (the local Mach number versus the normal coordinate y) on the 
symmetry plane at x = 0.04 (magenta line), x = 0.3 (red line) and x — 0.7 (blue line). The 
boundary 7 layer profile at the point (x,z) = (0.7,0.2) is shown by the green line, In all cases, 
there are more than 25 grid nodes in the boundary layer that provides sufficient accuracy for 
simulations of the near-wall How structures. 

In order to validate the code and grid refinements, we simulated the case of a — 20° without 
SDBD forcing and compared the numerical solution with that obtained with the How symmetry 
conditions. As shown in Fig. 33, the agreement in terms of the lift coefficient is good. 



Figure 33 C } (a) without SDBD forcing calculated using the flow symmetry conditions (red line) and 
for the total wing (purple triangle). 


2.7.1 Solutions for the wing-apex actuator 


Numerical simulations of the asymmetric SDBD forcing were conducted for the configuration 
shown in Fig. 34, where the forcing region is indicated by the red strip. I he SDBD induced 
source terms for the x- and y-momentum equations are approximated by the analytical relations, 
which arc similar to those given in Section 2.2. 

FA^t) = ^(t)A x ,g^x)f(y) (2.43) 

F, </• 1) = <h (x)f(v) 


Here < p(l) is a step-function: <p(t.) = 0 
distribution is 

<7 *•(*)= 1 


for /.<(), ^(/)=1 for />(). 

, , 0.2 

• x, < r < .r,. 


• r ~ - r i 

L, 


I he 


longitudinal 


(2.44) 


.t, = 2.5 cm, x., = .i\ f L,, where L, = 5 mm. I he distribution normal to the v\ ing surface is 
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/(?/) = a/i/exp 


(2.45) 




= 7547 7/* 


.(7y = 1() S /// 


!U = 0 - 1 , 


where the vertical size is approximately 0.4 mm, A rl , = 10 r> Njrn* and A ij{ , — 2 • 10 1 Njm . 


T he heat source is approximated as 

Q{d\yJ) ^ ^{t)A }1 (]{x)f{y) ^ 

where A ifT =10' W^m 3 . 


(2.46) 


Figure 35 shows 3D streamlines for symmetrical flow without SDBD forcing (a) and with the 
downstream forcing (b) for the angle of attack a = 20° . In this example, the heat source is zero. 
The streamlines are released from a small region near the wing apex. 



Figure 35 Streamlines and surface temperature at a = 20 . (a) SDBD is off (h) SDBD acts 
downstream on the bottom half plane, the heat source is zero. 
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It is seen that the vortex burst on the actuated wing surface is slightly delayed. On the non- 
actuatcd surface, the burst is shifted upstream significantly. This indicates that the SDBD effect 
is not local, and the wing-apex actuator affects the global How pattern in a nontrivial manner. 


Figure 36 illustrates the case when the SDBD actuator produces both the downstream 
momentum source and the heat source. It turned out that solutions were highly unsteady; i.e. the 
heat souree leads to significant increasing of the burst locus oscillations. 



Figure 36 Snapshots of 3D streamline pattern on the leeward wing surfaec. The SDBD actuator produces 
both the downstream momentum souree and the heat source on the bottom half plane. 

2.7.2 Solutions for the leading-edge actuator 

Consider the case when the SDBD region is located near the w ing leading edge as shown in Fig. 
37. The actuator ean bow out the wing symmetry plane or toward the wing symmetry plane as 
shown by the blue arrows. 



Figure 37 The leading-edge SDBD actuator blows out the wing symmetry plane (a) and toward the wing 
symmetry plane (b). 
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Figure 38 shows 3D streamlines over the leeward wing surfaee for the ease shown in Fig. 37a 
(the actuator blows away from the wing). The left and right plots correspond to different time 
instants. The SDBD region is schematically shown by the red dashed line. Sudden breakdowns 
of the streamline patterns indicate busts of the primary vortiees. 

Figure 39 illustrates the case when SDBD blows toward the wing symmetry plane (Fig. 37b). A 
qualitative behavior of the vortex How is similar to the previous ease The flow is highly 
unsteady. The vortex busts move baek and foree with appreciable amplitude. This leads to 
significant oscillations of the rolling moment 



Figure 38 3D streamlines over the leeward wing surface at different time instants. SDBD blows away 
from the wing (Fig. 37a), the actuator is schematically shown by the red dashed line 



Figure 39 3D streamlines over the leeward wing surface at different time instants. SDBD blows toward 
the wing (Fig. 37b), the actuator is schematically shown by the red dashed line. 
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2.7.3 Solutions for multi-element actuator 


Consider a multi-element SDBD schematically shown in Fig, 40. The red strips indicate the 
SDBD elements. The forcing is produced on one side of the wing leeward surface, where z > {). 
The coordinates x ui , x tfi are given in Table 3. The SDBD forcing is directed downstream. 



Figure 40 The multi-element actuator on the delta wing. The SDBD elements are shown by red lines. 


Table 3 Coordinates of SDBD elements 


/ = number of DBD strip 


x th , m 

1 

0.03 

0.045 

2 

0 15 

0.165 

3 

0.4 

0.415 

4 

0.6 

0.615 

5 

0.8 

0.815 


Figure 41 shows 3D streamline patterns over the leeward wing surface. The actuator is located 
on the low half plane. It is clearly seen that the SDBD forcing leads to appreciable delay of the 
vortex bust. Note that the vortex flow is highly unsteady similar to the eases considered in 
Sections 2.7.1 and 2.7.2 . The streamlines and the vortex bust loci oscillate that causes significant 
pulsations of the aerodynamic coefficients. 
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Figure 41 Streamline pattern over the leeward wing surface, the SDBD elements are located on the lower 
side (z>0) and blow downstream. 


3. Summary and Impact of Effort 

This project includes the two major components: modeling of the surface dielectric barrier 
discharge (SDE3D) and CFD studies of SDBD flow control capabilities. 

The first component has been focused on developments of self-consistent theoretical models and 
robust computational tools providing adequate from-first-principles simulations of heat and 
momentum sources induced by surface discharges. The physical model and numerical code were 
developed to predict SDBD process in atmospheric air. The model accounts for air 
photoionization by UV radiation from discharge zone, approximation of non local ionization 
source and universal boundary conditions on the electrode and dielectric surfaces. This allows us 
to simulate the total SDBD cycle including the discharge formation phase and the discharge 
relaxation phase. 

The numerical solutions showed that the discharge formation phase lasts few tens of 
nanoseconds and creates plasma. The discharge evolves as a streamer for positive electrode 
polarity and as a diffusive discharge for negative electrode polarity. The predicted values of the 
discharge length and the surface charge density agree well with the available experimental data 
both for positive and negative electrode polarity. This validates the developed physical model 
and computational code. 

Although the discharge formation phase gives negligible contribution to the momentum and heat 
sources relevant to flow control applications, it provides initial conditions for the relaxation 
phase. The latter begins when plasma is redistributed in such a way that it shields the external 
electric field and reduces the air ionization rate in the discharge region to very small value. The 
ion-electron and ion-ion recombination along with the ion drift motion are main processes in the 
relaxation phase, which lasts few microseconds. This phase effectively contributes to the 
momentum and heat sources required for flow control. The momentum-source distributions have 
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a complicated structure in space. For negative electrode polarity, the body force is an order of 
magnitude greater than in the case of positive applied voltage. 

Using these Findings and available experimental data we have developed a physics-based 
phenomenological model, which sheds light on processes governing the ignition and ceasing of 
microdischarges for the negative going half-cycle of applied alternating voltage. This model 
clarifies mechanisms leading to the body force saturation observed in experiments. The obtained 
analytical formulae capture dependencies of the body force and discharge length on the SDBD 
parameters. The model has been validated against the experimental data for SDBD actuators 
having dielectric layers of different materials and thicknesses. The phenomenological model led 
us to the assumption that the SDBD performance can be improved dramatically by increasing of 
the volumetric charge of negative ions. This can be done with the help of ‘conducting’ 
dielectrics. 

Using the body-force values predicted by the phenomenological model and numerical results for 
the spatial distributions of momentum and heat sources, we have obtained analytical 
approximations of the SDE^D-induced body force and heat source. The source terms were 
incorporated into Navier-Stokes solver. This computational tool has been used in the second 
component of the project dealing with SDBD applications to How control. 

To clarify aerodynamic effects produced by SDBD actuator, we have considered a relatively 
simple (low the laminar boundary layer on a flat plate in a free stream of low subsonic speed. 
Numerical solutions of 2D unsteady Navier-Stokes equations showed that the SDBD effect on 
the near-wall flow strongly depends on time-modulations of the applied voltage. If the discharge 
produces downstream momentum, then the actuator generates tangential near-wall jets. If the 
SDBD forcing is directed upstream, the actuator generates strong concentrated vortices. This 
flexibility and very short characteristic timescales make SDBD actuators attractive for control of 
the near-wall (lows associated with the boundary-layer separation, laminar-turbulent transition or 
a combination of both. 

CFD studies of the vortex flow' past a delta w ing at high angles of attack have been carried out to 
estimate feasibility of the vortex (low control using SDBD. Computations were performed for a 
delta wing with sharp leading edges of 60° sweep angle at the free-stream Maeh number 
A/ x = 1.5 , Reynolds number Rt\ = 2 x 10 f) and angles of attack 0 < cv < 30 . 

First series of computations was conducted without SDBD forcing. It was shown that our 
numerical solutions are consistent with available experimental and numerical results. Analysis of 
flow' patterns at various angles of attack indicated that the wing aerodynamic performance 
depends on the locus where the conical vortex structure breaks down. This suggested that the 
high-angle-of-attack flow control strategy should be focused on control of the vortex breakdown 
(or vortex burst) over the wing leeside. 

In the second series of computations, we have treated the wing-apex SDBD, leading-edge SDBD 
and multi-element SDBD actuators located on the wing leeside surface. The heat and momentum 
sources produced by these actuators were symmetric with respect to the wing centerline. It was 
found that the vortex-burst locus can be controlled by the aforementioned actuators. Namely, the 
breakdown point moves downstream with the SDBD acting in the downstream direction while it 
moves upstream with the SDBD acting upstream. However, the integral aerodynamic forces arc 
weakly affected. The actuator causes about 3% variations ot the lift and drag coefficients. 
Presumably the aerodynamic loads are weakly sensitive to the vortex burst locus for the delta 
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wing and frce-stream parameters considered herein. Direct numerical simulations of unsteady 
flow fields showed that the vortex breakdown evolves with time. Its unsteady behavior is 
sensitive to the SDBD forcing. Nevertheless, appreciable migrations of the burst locus produce 
small (less than 2%) effect on the lift coefficient. 

In the third series of computations, we considered the same SDBD actuators, which were located 
on one side of the wing surface. It was expected that an asymmetric SDBD forcing could 
generate appreciable rolling moment. However, numerical solutions of 3D Navier-Stokes 
equations indicated that the asymmetric forcing triggers strong oscillations of the vortex burst 
loci. This, in turn, causes pulsations of the rolling moment. Because of this detrimental effect it is 
not clear if SDBD actuators are robust for the vortex How control on supersonic delta wings at 
high angles of attack. Further numerical and experimental studies are needed to clarify this issue. 
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